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Abstract 

We investigate and compute spectral sums of the Wilson lattice Dirac operator for quenched SU(3) 
gauge theory. It is demonstrated that there exist sums which serve as order parameters for the confinement- 
deconfinement phase transition and get their main contribution from the IR end of the spectrum. They are 
approximately proportional to the Polyakov loop. In contrast to earlier studied spectral sums some of them 
are expected to have a well-defined continuum limit. 
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INTRODUCTION 



Confinement and chiral symmetry breaking are prominent features of strongly coupled gauge 
theories. If the gauge group contains a non-trivial center Z, then the traced Polyakov loop tii|2|] 

L x = tx c V x , V x = Y[U (t,x) (1) 

T = l 

serves as an order parameter for confinement in pure gauge theories or (supersymmetric) gauge 
theories with matter in the adjoint representation. The dynamics of L x near the phase transition 
point is effectively described by generalized Potts models yl, y|. Here we consider the space- 
independent expectation values (L x ) only and thus may replace L x by its spatial average 

L = 7tJ2 L *> V^Nt 1 . (2) 

V S x 

The expectation value (L) is zero in the center- symmetric confining phase and non-zero in the 
center-asymmetric deconfining phase. 

Chiral symmetry breaking, on the other hand, is related to an unusual distribution of the low 
lying eigenvalues of the Euclidean Dirac operator X> Q5Q. In the chirally broken low-temperature 
phase the typical distribution is dramatically different from that of the free Dirac operator since 
a typical level density p(A) for the eigenvalues per volume does not vanish for A — > 0. Indeed, 
according to the celebrated Banks-Casher relation flg], the mean density in the infrared is propor- 



tional to the quark condensate, 



(p(0)> = -i(0|^|0). (3) 

7T 



Which class of gauge field configurations gives rise to this unusual spectral behavior has not been 
fully clarified. It may be a liquid of instanton-type configuration QVD - Simulations of finite tempera- 
ture SU (3) gauge theory without dynamical quarks reveal a first order confinement-deconfinement 
phase transition at 260 MeV. At the same temperature the chiral condensate vanishes. This indi- 
cates that chiral symmetry breaking and confinement are most likely two sides of a coin ([8], for a 
review see e.g. D90). 

Although it is commonly believed that confinement and chiral symmetry breaking are deeply 
related, no analytical evidence of such a link existed up to a recent observation by Christof Gat- 
tringer lioll . His formula holds for lattice regulated gauge theories and is most simply stated for 
Dirac operators with nearest neighbor interactions. Here we consider fermions with ultra- local 



and 75-hermitean Wilson-Dirac operator 



d-l 



fi=0 



V\x) = (m + d)5 xy - - ((1 + 7 M )£Mz)<W-e M + (! - 7 m )^)<Wm 



(4) 



where U± ll (x) denotes the parallel transporter from site x to its neighboring site x ± such that 
U-^x + e^jU^x) = 1 holds true. Since we are interested in the finite temperature behavior we 
choose an asymmetric lattice with N t sites in the temporal direction and N s ^> N t sites in each of 
the d — 1 spatial directions. We impose periodic boundary conditions in all directions. The 



dim(D) = V x 2 [d/2] x N c , V = N t x 



(5) 



eigenvalues of the Dirac operator in a background field {U^(x)} are denoted by X p , The non-real 
ones occur in complex conjugated pairs since V is 75-hermitian. If N t and A^ s are both even, then 



A, 



2(d + m) — \ p is a further symmetry of the spectrum. 



x 



zUq(x) 

T 



Following [fiol . [llll we twist the gauge field 
configuration with a center element as fol- 
lows: all temporal link variables U (t, x) at 
a fixed time r are multiplied with an element 
z in the center Z of the gauge group. The 
twisted configuration is denoted by { Z U}. The 
Wilson loops Wc for all contractable loops C 
are invariant under this twisting whereas the 
Polyakov loops V x pick up the center element, 



W C ( Z U) = W C (U) and V X { Z U) = zV x {U). 



(6) 



The Dirac-eigenvalues for the twisted configuration are denoted by Z X P . The remarkable and simple 
identity in [10, 11] relates the traced Polyakov loop L to a particular spectral sum, 

1 \Z\ dim(X>) 



L = - £ z k 

k k=i p=i 



A, 



vJVi 



.l)^t 2 [d/2]-iy|2:|. 



(7) 



The first sum extends over the elements Zi,Z2,--- in the center Z containing the group identity 
e for which e X p = X p . The second sum contains the A^ t 'th power of all eigenvalues of the Dirac 
operator Zk V with twisted gauge fields { Zk U}. It is just the trace or ( Zk V) Nt , such that 



L = -J2z k tr( Zk V) 
K k 



Nt 



(8) 



We stress that the formula © holds whenever the gauge group admits a non-trivial center. In [fioll 



it was proved for SU(N C ) with center Z(iV c ) and k = |(— ) Nt dimV. In 111 ill the Dirac operator 
for staggered fermions and gauge group SU(3) was investigated and a formula similar to © was 
derived. Note that ® is not applicable to the gauge groups G 2 , F± and E 8 with trivial centers. 

For completeness we sketch the proof given in jtlOQ , slightly generalized to incorporate all 
gauge groups with non-trivial centers. The Wilson-Dirac operator contains hopping terms between 
nearest neighbors on the lattice. A hop from site x to its neighboring site x ± is accompanied 
by the factor — 1(1 =f 7 M )[/ At (x) and staying at x is accompanied by the factor m + d. Taking the 
£'th power of V, the single hops combine to chains of £ or less hops on the lattice. In particular the 
trace tr V e is described by loops with at most £ hops. Each loop C contributes a term proportional 
to the Wilson loop Wc- 

On an asymmetric lattice with N t < N s all loops with length < N t are contractable and since 
the corresponding Wilson loops Wc do not change under twisting one concludes 

tr z X? £ = triy for £ < N t . (9) 

For any matrix group with non-trivial Z the center elements sum to zero, J2 %k = 0, such that 

5> fc tr( 2 ^)=tr(:r/)E^ = for £ < N t . (10) 

k k 

For £ = N t only the Polyakov loops winding once around the periodic time direction are not 
contractable. Under a twist by {U} — > { Z U} they are multiplied by z, see ©. With J2k %kZk = \Z\ 
we end up with the result © which generalizes Gattringer formula to arbitrary gauge groups with 
non-trivial center. What happens for £ > N t in (flOl) will be discussed below. 



In ||1 111 the average shift of the eigenvalues when one twists the configurations has been cal- 
culated. It was observed that above T c the shift is greater than below T c and that the eigenvalues 
in the infrared are more shifted than those in the ultraviolet. But the low lying eigenvalues are 
relatively suppressed in the sum © such that the main contribution comes from large eigenvalues. 
Indeed, if one considers the partial sums 

X„ = -£^E^f, n<dim(V), (11) 

K k p=l 

where the eigenvalues are ordered according to their absolute values, then on a 4 3 x 3-lattice 70% 
of all eigenvalues must be included in (fTTT) to obtain a reasonable approximation to the traced 



Polyakov loop llllh . Actually, if one includes fewer eigenvalues then the partial sums have a phase 



shift of 7r relative to the traced Polyakov loop. For large N t the contribution from the ultraviolet 
part of the spectrum dominates the sum ©. Thus it is difficult to see how the nice lattice result © 
could be of any relevance for continuum physics. 

The paper is organized as follows: In the next section we introduce flat connections with zero 
curvature but non-trivial Polyakov loops. The corresponding eigenvalues of the Wilson-Dirac op- 
erator are determined and spectral sums with support in the infrared of the spectrum are defined 
and computed. The results are useful since they are in qualitative agreement with the correspond- 
ing results of Monte-Carlo simulations. In section 3 we recall the construction of the real order 
parameter L mt related to the Polyakov loop [4]. Its Monte-Carlo averages are compared with 
the averages of the partial sums (fTTI) . Our results for Wilson-Dirac fermions are in qualitative 
agreement with the corresponding results for staggered fermions in illlh . In section 4 we discuss 
spectral sums for inverse powers of the eigenvalues. Their Monte-Carlo averages are proportional 
to (L rot ) such that they are useful order parameters for the center symmetry. We show that these 
order parameters are supported by the eigenvalues from the infrared end of the spectrum. Section 
5 contains similar results for exponential spectral sums. Again we find a linear or quadratic rela- 
tion between their Monte-Carlo averages and (L TOt ). It suffices to include only a small number of 
infrared eigenvalues in these sums to obtain efficient order parameters. We hope that the simple 
relations between the infrared- supported spectral sums considered here and the expectation value 
(L rot ) are of use in the continuum limit. 



FLAT CONNECTIONS 



We checked our numerical algorithms against the analytical results for curvature-free gauge 
field configurations with non-trivial Polyakov loop. For these simple configurations the spatial 
link variables are trivial and the temporal link variables are space-independent, 

Ui(x) = l and U q (x) = U (t), x = (t,x). (12) 

The Wilson loops Wc of all contractable C are trivial which shows that these configurations are 
curvature-free. We call them flat connections. With the gauge transformation 

Q(r) = V-\ V T = U (t - 1)U (t - 2) • • • U (2)U (1) (13) 
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all link-variables of a flat connection are transformed into the group-identity. But the transformed 
fermion fields are not periodic in time anymore, 

V>(r + jVt,x) =V- 1 i(){t,x), where V = V Nt+1 (14) 

is just the constant Polyakov loop. Since the transformed Dirac operator is the free operator, its 
eigenfunctions are plane waves, 

xl>(x) = e^Vo- (15) 
These are eigenmodes of the free Wilson-Dirac operator with eigenvalues {A p } = {A^ }, where 

\p=m±i\p\ + — , with ^ = 2sin^, p^ = sinp M . (16) 
They are periodic in the space directions provided the spatial momenta are from 

Pi e — rii with m e Zjv s . (17) 

Denoting the eigenvalues of the Polyakov loop by e 2nllpl , . . . , e 2m<fiNc , the periodicity conditions 
(TT41) imply 

Po = — (n - <fj), n eZ Nt , j = l,...,N c . (18) 

Thus the eigenvalues of the Wilson-Dirac operator with a flat connection are given in (fT6l) . with 
quantized momenta (TT71 ) and (TT81 ). For each momentum there exist 2^ 2 ^ 1 eigenvalues A+ and 
2[ d / 2 M complex conjugated eigenvalues A~. 

Next we twist the flat connections with a center-element, for SU(N C ) with 



z k 



e 2nik / N "l, l<k<N c . (19) 



The spatial components of the momenta are still given by (fTTT) . but their temporal component is 
shifted by an amount proportional to k, 

2tt 



Po(z*) 6 {-^ (n - ^ - fe/iV c )| , 1 < j, fc < jV c . 



(20) 



In the following we consider flat S'f/(3)-connections with Polyakov loops 

/ e 2ni8 q \ 

V(0) =010 = 
V e~ 2md / 



L = 1 + 2cos(27r#). (21) 



For these fields the temporal component of the momentum takes values from 

Po(**)e{|^(no-j0-V3)}, 3 e {-1,0,1}, k G {0, 1, 2}. 



(22) 



We have calculated the spectral sums 

i dimZ> i 

s w = -£** £ = -£** tr (23) 

K it p=i K fc 

for vanishing mass. For flat connections the sums with powers £ between iVt and 2iV t are strictly 
proportional to the traced Polyakov loop, = CgL(9). Gattringers result implies Cjsr t = 1. The 
next two coefficients are related to the number of loops of length N t + 1 and iV t + 2 winding once 
around the periodic time direction. One finds 

C Nt+1 = d(N t + l) and C Nt+2 = ^(N t + 2)(N t + l) + ^—±(N t (N t + l)-2). (24) 

More generally, the relation ^k z i = for £ ^ 3Z + 1 implies that the spectral sums (|23T) are 
linear combinations of the traces tr'P 3n+1 (6 l ) for sufficiently small values of \3n + 1|, 



£ , 

n: |3n+l|7V t <£ 



d n) trP 



3n+l 



(9). 



(25) 



In Fig. [Qwe depicted the sums on a 4 x 12 3 lattice, divided by the traced Polyakov loop and 
normalized to one for 9 = for the flat connections and the powers I — N t , 3N t and 3.QN t . Note 
that the power £ in (|23l) need not be an integer. 



1.00 



£ = N t 




e = 3-N t 




£ = 3.6- N t 



9 



0.95 

1 
FIG. 1: Spectral sums divided by the traced Polyakov loop as functions of 9 for different values of £. 

We have argued that the sum £W must be a linear combination of tr V and tr V~ 2 for £ between 
2N t and AN t . Actually, up to £ « 3iV t the sum is well approximated by a multiple of tr V. This 
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is explained by the fact that for a given £ there are much more fat loops winding once around the 
periodic time direction and contributing with tr V than there are thin long loops winding many 
times around and contributing with tr V~ 2 , tr V 4 , trV~ 5 , . . .. We shall see that similar results 
apply to the expectation values of in Monte-Carlo generated ensembles of gauge fields. 

Since the eigenvalues in the infrared are mostly affected by the twisting [ljj] we could as well 
choose a spectral sum for which the ultraviolet end of the spectrum is suppressed. Since with 
£ < 3N t is almost proportional to the traced Polyakov loop there exist many such spectral sums. 
They define order parameters for the center symmetry and may possess a well-defined continuum 
limit. For example, the exponential sums 

dim© 

(26) 



5>* E e 

k p=l 



-e-x p ( z kU) 



are all proportional to the traced Polyakov loop for a factor £ in the exponent between 0.1 and 2. 
Below we displayed exponential sums for the flat connections on a 4 x 12 3 -lattice and various £ 
between 0.1 and 2. Again we divided by the traced Polyakov loop L(9) and normalized the result 
to unity for 9 = 0. 
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FIG. 2: Spectral sums divided by the traced Polyakov loop as functions of 6 for different values of I. 

Later when we use Monte-Carlo generated configurations to calculate the expectation values of L 
and £w we shall choose £ = 1. For this choice the mean exponential sum will be proportional to 
the mean L. Later we shall argue why this is the case. 
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DISTRIBUTION OF DIRAC EIGENVALUES FOR SU(3) 



We have undertaken extended numerical studies of the eigenvalue distributions and various 
spectral sums for the Wilson-Dirac operator in SU (3) lattice gauge theory. First we summarize 
our results on the partial traces 

1 



■£>£*Aj> n<dim(V), 



< |A p +i| 



(27) 



k p=l 

For n = dim(V) one sums over all eigenvalues of the Dirac-operator and obtains the traces 
considered in (T23T ). For i = N t one finds the partial sums S n in (fTTI) . These have been extensively 
studied for staggered fermions in llllll . According to the result © the object £ dim x> is just the 
traced Polyakov loop. 

We did simulations on lattices with sizes up to 8 3 x 4. Here we report on the results obtained 
on a 4 3 x 3 lattice with critical coupling /3 crit ~ 5.49, determined with the histogram method based 
on 40 000 configurations. The dependence of the two order parameters \L\ and U ot (see below) on 
the Wilson coupling (3 has been calculated for 35 different (3 and is depicted in Fig. |3] For each (3 

1.4 




FIG. 3: Dependence of the mean modulus of L and the center-transformed and rotated L (see text) on the 
Wilson coupling (3 on a 4 3 x 3 lattice. The critical coupling is f3 CI n = 5.49 



between 4 000 and 20 000 independent configuration have been generated and analyzed. For our 
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relatively small lattices the two order parameters change gradually from the symmetric confined 
to the broken deconfined phase. Table U contains the order parameters for 11 Wilson couplings. 
For every independent configuration we calculated the dim(D) = 2304 eigenvalues of the Wilson- 



p 


5.200 5.330 5.440 5.475 


5.505 5.530 5.560 5.615 5.725 5.885 6.000 


(\L\) 
(L rot ) 


0.1654 0.1975 0.3050 0.3980 
0.0318 0.0615 0.1859 0.3013 


0.5049 0.5939 0.6865 0.7832 0.9007 1.0013 1.0631 
0.4296 0.5363 0.6452 0.7513 0.8770 0.9797 1.0444 



TABLE I: Dependence of the order parameters \L\ and L rot on the Wilson coupling (3. 



Dirac operator. Then we averaged the absolute values of the partial traces S n for every (3 in table 
HI In Fig. @]the ratios 

<28) 

are plotted for these (3 as function of the percentage of eigenvalues considered in the partial traces. 



1.4 




% of lowest eigenmodes 



FIG. 4: Modulus of the eigenvalue sums starting from the lowest eigenmodes on a 4 3 x 3-lattice near the 
phase transition. The distinct graphs are labelled with the Wilson coupling 0. 



To retain information on the phase of the partial traces and Polyakov loop we used the invariant 
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order parameter constructed in |4j]. Recall that the domain for the traced Polyakov loop is the 
triangle shown in Fig. [5] The three elements in the center of £77(3) correspond to the corners of 
the triangle. We divide the domain into the six distinct parts in Fig. |5J The light-shaded region 




FIG. 5: Fundamental domain T of L obtained by identifying Z(3) copies according to the depicted arrows. 



represents the preferred location of the traced Polyakov loop L in the deconfined (ferromagnetic) 
phase, whereas the dark-shaded region corresponds to the hypothetical anti-center ferromagnetic 
phase 1120. In the deconfined phase L points in the direction of a center element whereas it points 
in the opposite direction in the anti-center phase. To eliminate the superfluous center- symmetry we 
identify the regions as indicated by the arrows in Fig. \5\ This way we end up with & fundamental 
domain T for the center-symmetry along the real axis. Every L is mapped into T by a center 
transformation. To finally obtain a real observable we rotate the transformed L inside T onto the 
real axis. The result is the variable L TOt whose sign clearly distinguishes between the different 
phases. L TOt is negative in the anti-center phase, positive in the deconfined phase and zero in 
the confined symmetric phase. The object L rot is a useful order parameter for the confinement- 
deconfinement phase transition in gluodynamics 04]. 

We performed the same construction with the partial sums E n and calculated the ratios for the 
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corresponding Monte-Carlo averages 



Rl 



(L rot ) 



(29) 



for every /3 in table U as a function of the percentage of eigenvalues considered in E n . In Fig. @]and 
1.0 
0.8 




30 40 50 60 70 
% of lowest eigenmodes 



100 



FIG. 6: Eigenvalue sums rotated to the fundamental domain starting from the lowest eigenmodes on a 
4 3 x 3-lattice near the phase transition. The distinct graphs are labelled with the Wilson coupling. 

|6] we observe a universal behavior in the deconfined phase with modulus of the traced Polyakov 
loop larger than approximately 0.4. If we include less than 90% of the eigenvalues, then the partial 
sums E n have a phase shift of n in comparison with E = E dim £). The last dip in Fig|4]is due to this 
phase shift and indicates the transition through zero that occurs when E n changes sign. The same 
shift and dip has been reported for staggered fermions on a 6 3 x 4 lattice in 111 ill . For staggered 
fermions E n and E are in phase for n > 0.65 ■ dim(P). For Wilson-Dirac fermions this happens 
only for n > 0.9 ■ dim(P). 

Finite spatial size scaling of partial sums: We fixed the coupling at (5 = 5.5 and simulated 
in the deconfined phase on iV s 3 x 2-lattices with varying spatial sizes iV s 6 {3,4,5}. For this 
coupling the systems are deep in the broken phase and we can study finite size effects on the 
spectral sums. The results for E™* are depicted in Fig. [71 We observe that to a high precision 
E™* is approximately independent of the spatial volume. The curves for iV s = 4 and 5 are not 
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Tot\ 




30 40 50 60 70 
% of lowest eigenmodes 



100 



FIG. 7: Rotated eigenvalue sums starting from the lowest X p on a iV s 3 x 2-lattice in the broken phase. 

distinguishable in the plot and as expected J2k Zfctr(**P *) scales with the spatial volume of the 
system. An increase of N s affects the spectra for the untwisted and twisted configurations alike - 
they only become denser with increasing spatial volume. On the other hand, comparing Fig. [6]and 
Fig. |7J it is evident that the graph of E™* depends very much on the temporal extent of the lattice. 

Partial traces E^; The truncated eigenvalue sums (|27T) with different powers £ of the eigen- 
values show an universal behavior that is nearly independent of the lattice size. The main reason 
for this universality and in particular the sign of is found in the response of the low-lying 
eigenvalues to twisting the gauge field. It has been observed that for non-periodic boundary condi- 
tions (which are gauge-equivalent to twisting the gauge field) the low lying eigenvalues are on the 
average further away from the origin as compared to periodic boundary conditions (or untwisted 



gauge fields) [13 



1J, 



15 



This statement is very clear for massless staggered fermions with 
eigenvalues on the imaginary axis. For example, the partial traces 

\\p\ < \\p+i\ 



EWoc^Ap + ^E^ + ^E^ 
p=i p=i p=i 



V 



(30) 



with n <C dim V and the traced Polyakov loop have opposite phases. This is explained as follows: 
all sums in ( |30l ) are positive and on the average the last two sums are equal. With z + z = — 1 the 
last two terms add up to — X) z \ p . Since the low lying eigenvalues for the twisted field are further 
away from the origin as for the untwisted field, the spectral sums (l30l) are negative for small n. 
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TRACES OF PROPAGATORS 



To suppress the contributions of large eigenvalues we introduce spectral sums with negative 
exponents L Similar to the Polyakov loop these sums serve as order parameters for the center 
symmetry. In particular the spectral sums 

S (-D = IVtr(4^) and £(- 2 ) = -Vtr(^V) (3D 

are of interest, since they relate to the Green functions of V and V 2 , objects which enter the 
discussion of the celebrated Banks-Casher relation. Contrary to the ultraviolet-dominated sums 
with positive i are the sums with negative t dominated by the eigenvalues in the infrared. The 
operators { Z V, z G Z} have similar spectra and we may expect that kTS~^ has a well-behaved 
continuum limit. Here we consider the partial traces 

Since the Wilson-Dirac operator with fiat connection possesses zero-modes we added a small 
mass m = 0.1 to the denominators in (|32l) . In Fig. [8] the partial sums S^ -1 - 1 on a 4 3 x 3 lattice 
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1 1 1 



25 50 75 100 



% of lowest eigenvalues 
FIG. 8: The partial spectral sums ^ for the inverse power and flat connections. 

are plotted. It is seen that for flat connections the E^ -1 ) for small n are excellent indicators for 
the traced Polyakov loop. Thus it is tempting to propose S^ - ^, S^~ 2 ) with n <^ dimD as order 
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parameters for the center symmetry. To test this proposal we calculated the partial sums (|32|) . 
transformed to the fundamental domain and rotated to the real axis, for Monte-Carlo generated 
configurations on a 4 3 x 3 lattice for various values of (3. The results in Fig. [9] are qualitatively 
similar to those for the fiat connections. Taking into account 10% of the eigenvalues in the IR 
already yields the asymptotic values £(~ 1 )' rot and E( _2 )> rot . 

Or 

xlO" 3 



-,(-l),rat 




10 20 30 40 50 60 70 80 90 100 



xlO" 



v (-2),rot 




10 20 30 40 50 60 70 80 90 100 
% of lowest eigenmodes 



FIG. 9: The expectation values of the partial spectral sums S„ ^ and En ^ rotated to the fundamental 
domain starting from the lowest eigenvalue on a 4 3 x 3 lattice. The graphs are labelled with /?. 
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To find an approximate relation between E^ 1 ) and the traced Polyakov loop we applied the 
hopping-parameter expansion. To that end one expands the inverse of the Wilson-Dirac operator 
V = (m + d) 1 — V in powers of V, 

1 ^ 1 



V 



-i 



y 

m + d ^ (m + dy 



■[(m + d)l-VY 



(33) 



Inserting this Neumann series into E^ ^ in (1311) and keeping the leading term only yields 



S (-D 



■1) 



E« + . . . 



(8, 



-L. 



(34) 



To check whether the expectation values of E( _1 ^ ,rot and L rot are indeed proportional to each 
other we have calculated these values for Monte-Carlo ensembles corresponding to the 1 1 Wilson 
couplings in table[D The results in Fig. [101 clearly demonstrate that there is such a linear relation. 



xlO" 3 



(£(-!W) 




FIG. 10: The expectation values of £( -1 ) ,rot as functions of (L rot ) on a 4 3 x 3 lattice. 



A linear fit yields 



( S (-i),™t) = _ .00545 ■ (L rot ) - 4.379 ■ KT (rmse = 2.978 ■ 10 



(35) 

For massless fermions on a 4 3 x 3 lattice the crude approximation (|28l) leads to a slope —0.003906. 
This is not far off the slope —0.00545 extracted from the Monte-Carlo data. 
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We have repeated our calculations for the spectral sum £(~ 2 ) rot . The corresponding results for 
the expectation values in Fig. \TT\ show again a linear relation between the expectation values of 
this spectral sum and the traced Polyakov loop. 



xlO" 3 



/ E (-2),rotN 




FIG. 11: The expectation values of £( 2 '' rot as functions of (L rot ) on a 4 3 x 3 lattice. 



This time a linear fit yields (S(~ 2 )' rot ) = -0.00582 • (L rot ) - 8.035 • 10~ 5 . 



EXPONENTIAL SPECTRAL SUMS 



After the convincing results for sums of inverse powers of the eigenvalues we analyze the partial 
exponential spectral sums 



n 

k p=l 
n 



Qn = - E Zk E ' 



p=l 



£ = SdimV = ~ E Z k tr eX P (- 2fe ^) 
K k 

K. * — ' V 



(36) 
(37) 



In particular the last expression is used in a heat kernel regularization of the fermionic determinant. 
kQ has a well-behaved continuum limit if we enclose the system in a box with finite volume. We 
computed the partial sums Q n for the flat connections and various values of the traced Polyakov 
loop. In Fig. [12] we plotted those sums for which 10% or less of the low lying eigenvalues have 
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FIG. 12: The partial Gaussian sums Q n for fiat connections with different L. 

been included. Similarly as for the sums of negative powers of the eigenvalues we conjecture that 
the Gaussian sums Q n are good candidates for an order parameter in the infrared. 

The expectation values of the partial sums £* ot and Q™* for Monte-Carlo generated configura- 
tions at four Wilson couplings are plotted in Fig. [13] and Fig. [T4] As expected from our studies 




10 20 30 40 50 60 70 80 90 100 



% of lowest eigenmodes 

FIG. 13: Mean exponential sums ££ ot on a 4 3 x 3-lattice near /? cr it- The graphs are labelled with (3. 

of fiat connections, the Gaussian sums are perfect order parameters for the center symmetry. They 
are superior to the other spectral sums considered in this paper, since their support is even further 
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at the infrared end of the spectrum. Fig. [15] shows the expectation values (G™*) with only 4.5% or 
less of the infrared-modes included. The result is again in qualitative agreement with that for fiat 

0.0 -. — 
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FIG. 14: Mean Gaussian sums Q T ° l on a 4 3 x 3-lattice near /3 crit . The graphs are labelled with (3. 
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FIG. 15: Zooming into Gaussian sums £7£ ot on a 4 3 x 3-lattice near the phase transition. 



connections in Fig. [121 although in the Monte-Carlo data the dips are washed out. 
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The Monte-Carlo results for the expectation values (£ rot ) and (L rot ) with Wilson couplings in 
table U are depicted in Fig. [TH The quality of the linear fit 



(S TOt ) = -0.00408 • (L rot ) + 2.346 • 10~ 5 (raise = 1.82 • 10~ 5 ), 



is as good as for the spectral sum E^ 1 \ 

1 

xlO" 3 




(38) 



FIG. 16: The expectation value of £ rot as function of (L rot ) on a 4 3 x 3 lattice. 



To estimate the slope and in particular its dependence on the lattice size we expand the exponentials 
in £ n which results in 



(39) 



Since is proportional to the traced Polyakov loop for i < 3iVt we conclude that £ = £dimT> 
should be proportional to L. We can estimate the proportionality factor as follows: in the Wilson 
loop expansion of trV^ Nt+p ^ only loops winding around the periodic time direction contribute. 
If we neglect fat loops and only count the straight loops winding once around the periodic time 
direction, then there are 



(m + d) 1 



P 



(40) 
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such loops contributing. Actually, for p > N t there are loops winding several times around the 
time direction. But these have relatively small entropy and do not contribute much. Hence, with 
(1391 ) we arrive at the estimate 

S - (-U* £ <^> H- W . ") . L = <=gVo~ ... (41, 

In 4 dimensions and for m = we obtain the approximate linear relation 

N t \£ « (-1) M -0.0183 -L. (42) 

For the linear fit ([38]) to the MC-data the slope is 3! • 0.00408 = 0.0245 instead of 0.0183. 

The Monte-Carlo results for the order parameters (G rot ) and (L rot ) with Wilson couplings from 
table U are shown in Fig. [FT] In this case the functional dependence is more accurately described 



(£ rot ) 




0.2 0.4 0.6 0.8 1 



(L rot ) 

FIG. 17: The expectation value of Q mt as function of (L rot ) on a 4 3 x 3 lattice. 

by a quadratic function, 

(g mt ) = -0.000571 • (L rot ) 2 - 0.00156 • (L rot ) + 1.061 • 10~ 5 (rmse = 1.453 • 10" 5 ), (43) 

and this relation is very precise. Since in addition (Q^) ~ (G TOt ) already for small n we can 
reconstruct the order parameter (L rot ) from the low lying eigenvalues of the Wilson-Dirac operator. 
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Scaling with N t : On page [12] we discussed the finite (spatial) size scaling of the MC expec- 
tation values (S^ ot ). We showed that they converge rapidly to their infinite- iV s limit, see Fig. [7] 
Here we study how the Gaussian sums Q n depend on the temporal extend of the lattice. To that 
end we performed simulations on larger lattices with fixed iV s = 6, variable N t = 2, 3, 4, 5 and 
Wilson coupling (3 = 6.5. We calculated the ratios 



RL 



(L rot ) 



(44) 



where we multiplied with the extensive factor k in © since in the partial sums 

Gn 



k p=l 



\\\ <■ l-V+il- 



(45) 



only a tiny fraction of the 5184 to 12 960 eigenvalues have been included. The order parameter 
(L rot ) for the lattices with N t = 2, 3, 4, 5 is 1.9474, 1.40194, 0.932245, 0.523142. In Fig. [18]we 
plotted the ratios R™* for n from 1 up to 100. Note that on the 6 3 x 5-lattice n = 100 means less 

9 



R 



rot 




FIG. 18: The rations i^ ot as function of the number n of IR-eigenvalues included. 100 eigenvalues corre- 
sponds to approximately 1% of all eigenvalues. 



than 0.8% of all 12 960 eigenvalues. This figure very much supports our earlier statements about 
the quality of the order parameters (G%*) or (G^). 
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CONCLUSIONS 



In this paper we studied spectral sums of the type 

in in 

£n(/) = -E^E/r%) and 4(/) = ~E^E/(r A /)- ( 46 ) 



K k p=l K k p=l 



where { Zfc A p } is the set eigenvalues of the Wilson-Dirac operator with twisted gauge field. Sum- 
ming over all dim(£>) eigenvalues the sums over p become traces such that 

s(f) = - E tr / ( ZkV ) and ^(/) = - E tr / • ( 47 ) 

K k K k 



For /(A) = A^* one finds the spectral sum £ which reproduces the traced Polyakov loop [10]. 
Unfortunately this lattice-result is probably of no use in the continuum limit. Thus we have used 
functions /(A) which vanish for large (absolute) values of A. The corresponding sums are order 
parameters which get their main contribution from the infrared end of the spectrum. Of all spectral 
sums considered here, the Gaussian sums Q n in (|37l) define the most efficient order parameters. 
Besides the Q n the spectral sums of inverse powers of eigenvalues are quite useful as well. This 
observation may be of interest since these sums relate to the Banks-Casher relation. 

It remains to investigate the continuum limits of the spectral sums considered in this paper. The 
properly normalized Q n should have a well-behaved continuum limit. With regard to the conjec- 
tured relation between confinement and chiral symmetry breaking it would be more interesting 
to see whether the suitably normalized sums E^ 1 ) or/and S^ 2 ^ can be defined in the continuum 
theory. Clearly, the answer to this interesting question depends on the dimension of spacetime. 

Acknowledgments: We thank Georg Bergner, Falk Bruckmann, Christof Gattringer, Tobias 
Kastner and Sebastian Uhlmann for interesting discussions. This project has been supported by 
the DFG, grant Wi 777/8-2. 



[1] A. M. Polyakov, Phys. Lett. B72, 477 (1978). 

[2] L. Susskind, Phys. Rev. D20, 2610 (1979). 

[3] L. G. Yaffe and B. Svetitsky, Phys. Rev. D26, 963 (1982). 

[4] C. Wozar, T. Kaestner, A. Wipf, T. Heinzl, and B. Pozsgay, Phys. Rev. D74, 114501 (2006), hep- 
lat/0605012. 

[5] H. Leutwyler and A. Smilga, Phys. Rev. D46, 5607 (1992). 

23 



[6] T. Banks and A. Casher, Nucl. Phys. B169, 103 (1980). 

[7] T. Schafer and E. V. Shuryak, Rev. Mod. Phys. 70, 323 (1998), hep-ph/96 10451. 

[8] J. B. Kogut et al., Phys. Rev. Lett. 50, 393 (1983). 

[9] F. Karsch, Lect. Notes Phys. 583, 209 (2002), hep-lat/0106019. 

[10] C. Gattringer, Phys. Rev. Lett. 97, 032003 (2006), hep-lat/0605018. 

[11] F. Bruckmann, C. Gattringer, and C. Hagen (2006), hep-lat/0612020. 

[12] A. Wipf, T. Kaestner, C. Wozar, and T. Heinzl, SIGMA 3, 006 (2007), hep-lat/06 10043. 

[13] S. Chandrasekharan and N. H. Christ, Nucl. Phys. Proc. Suppl. 47, 527 (1996), hep-lat/9509095. 

[14] S. Chandrasekharan and S.-z. Huang, Phys. Rev. D53, 5100 (1996), hep-ph/95 12323. 

[15] M. A. Stephanov, Phys. Lett. B375, 249 (1996), hep-lat/9601001. 

[16] C. Gattringer, P. E. L. Rakow, A. Schafer, and W. Soldner, Phys. Rev. D66, 054502 (2002), hep- 
lat/0202009. 



24 



